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We present an improved numerical relativity (NR) calibration of the new effective-one-body (EOB) 
model for coalescing non precessing spinning black hole binaries recently introduced by Damour and 
Nagar [Physical Review D 90 , 044018 (2014)]. We do so by comparing the EOB predictions to both 
the phasing and the energetics provided by two independent sets of NR data covering mass ratios 
1 < q < 9.989 and dimensionless spin range —0.95 < X < +0.994. One set of data is a subset of 
the Simulating eXtreme Spacetimes (SXS) catalog of public waveforms; the other set consists of 
new simulations obtained with the Llama code plus Cauchy Characteristic Evolution. We present 
the first systematic computation of the gauge-invariant relation between the binding energy and 
the total angular momentum, Eb(j), for a large sample of, spin-aligned, SXS and Llama data. The 
dynamics of the EOB model presented here involves only two free functional parameters, one (og (u)) 
entering the non spinning sector, as a 5PN effective correction to the interaction potential, and one 
(c3(5i, 02 , i')) in the spinning sector, as an effective next-to-next-to-next-to-leading order correction 
to the spin-orbit coupling. These parameters are determined (together with a third functional 
parameter AInqc (x) entering the waveform) by comparing the EOB phasing with the SXS phasing, 
the consistency of the energetics being checked afterwards. The quality of the analytical model 
for gravitational wave data analysis purposes is assessed by computing the EOB/NR faithfulness. 

Over the NR data sample and when varying the total mass between 20 and 200 Mq the EOB/NR 
unfaithfulness (integrated over the NR frequency range) is found to vary between 99.493% and 
99.984% with a median value of 99.944%. 

PACS numbers: 04.30.Db, 04.25.Nx, 95.30.Sf, 97.60.Lf 


I. INTRODUCTION 

The purpose of this paper is to present an improved 
version of the new spinning effective-one-body (EOB) 
model introduced in Ref. /T] , hereafter Paper I. [For ear¬ 
lier spinning EOB models see nm.] When the mass of 
the binary varies between 20 and 200 Mq , our improved 
model yields, when compared to a large sample of numer¬ 
ical relativity (NR) waveform data, maximal unfaithful¬ 
nesses (F = 1 — F, integrated over the NR frequency 
range) ranging between 0.00016 and 0.00507, with me¬ 
dian value equal to 0.00056. The impatient data-analyst 
reader will find this information in Figs. [9] and [2l] and 
the text around. The structural elements of our improved 
EOB model behind such good faithfulnesses are: (i) the 
incorporation of the full 4PN-accurate analytical knowl¬ 
edge of the EOB radial interaction; and (ii) the use of 
a recently proposed improved ringdown description |14| . 
The use of the 4PN information obliged us to update 
the calibration of both the effective 5PN-coefficient a% 
entering the nonspinning sector, and the effective next- 
to-next-to-next-to-leading order spin-orbit coupling coef¬ 
ficient C3 of Paper I. 

Here, we restrict attention to the nonprecessing case 
where the spins are either aligned or antialigned with re¬ 
spect to the orbital angular momentum. The EOB/NR 
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comparisons that are used to gain new information from 
NR data so as to complete the EOB model are of two 
kinds: (i) on the one hand, we compare the gravita¬ 
tional wave (GW) phasing of the EOB model with 40 
state-of-the-art (publicly available [TS]) NR waveforms 
produced by the Caltech-Cornell-CITA Simulating eX¬ 
treme Spacetimes (SXS) collaboration with the Spectral 
Einstein Code (SpEC) code [T6U25] : eleven configurations 
in this sample involve nonspinning binaries, while in the 
remaining 29 at least one black hole is spinning; (ii) on 
the other hand we compare the EOB and NR energetics 
through the gauge-invariant relation between energy and 
angular momentum. To do so, we employ, in addition to 
the SXS datasets, ten newly performed simulations ob¬ 
tained with the Llama code, as a follow up of our previous 
work [25] , 


The paper is organized as follows: in Sec. |TT]we briefly 
review the origin of the NR data used in this paper, which 
were obtained with very different codes. In Secs |III||TV| we 
focus on the improved calibration of the nonspinning sec¬ 
tor, its phasing and energetics, notably showing the good 
agreement between energetics when J~ r ^ = 0 for the mass 
ratio range 1 < q < 9.989. In Secs. |Vl|VT| we calibrate the 
effective spin-orbit parameter C3 and we assess the quality 
of this new spinning EOB model by computing both the 


EOB/NR unfaithfulness and the energetics. In Sec. VII 


we present a preliminary comparison of our EOB model 
with that of Alessandra Buonanno’s group Ej. Finally, 
Sec. |VIII| presents concluding remarks, future prospects 
as well as a histogram summarizing the unfaithfulness 
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calculations. 

We use geometrized units with G = c = 1 and the fol¬ 
lowing notation: M = TO1+TO2, /r = TO1TO2/M, v = fi/M 
with the convention that the mass ratio q = mi/m 2 > 1. 

II. NUMERICAL RELATIVITY DATA 
A. Data from the SXS public catalog 

We use a sample of 40 simulations from the SXS cata¬ 
log. We consider either nonspinning or spin-aligned con¬ 
figurations, with mass ratios in the range 1 < q < 9.989 
and dimensionless spin magnitude —0.95 < \ < +0.994. 
In several cases, just one of the two black holes is spin¬ 
ning. The complete information about the data we used 
is listed in Table [I] For each run, when it is possible, 
we add an error bar on phasing computed at the NR 
merger. Here and below the instant of “merger” is de¬ 
fined as the time of the location of the peak of the mod¬ 
ulus of the l = m = 2 waveform. This error bar was esti¬ 
mated by computing the value at (the highest-resolution- 
) merger of the phase difference between the highest (H) 
and second-highest (SH) resolutions available in the cat¬ 
alog, 5 NR (j)(t) = dP(t) — </> SH (f). This difference is mono- 
tonically varying N up to merger and, by definition, taken 
to be zero at the start of the two simulations. When only 
one resolution is available we do not indicate any phasing 
uncertainty in Table [T] In such cases one can check the 
consistency by looking at datasets with neighbouring pa¬ 
rameters. In SXS simulations, the data are extracted at 
finite radius and then extrapolated to future null infinity 
(e.g. [2Zi[28]). This is done by means of a fit using a poly¬ 
nomial in 1/r, where r is the extraction radius. The user 
of the catalog is free to choose between different orders of 
extrapolation (N = 2 ,3,4) depending on the application, 
with the warning that the “best” extrapolation order to 
use depends on the simulation and that higher extrap¬ 
olation orders tend to do better during the inspiral and 
worse during ringdown. We found experimentally that 
extrapolation order N = 3 is reliable for all datasets ex¬ 
cept for SXS:BBH:0002, where we used N = 2 to reduce 
unphysical oscillations during the late inspiral. In addi¬ 
tion to the data of the catalog, we also use a ~ 14 orbits 
q = 1 , x = 0 waveform [J35], that was used to calibrate 
earlier EOB models |5J [25]. 


B. Data from the Llama code 

We use a sample of 10 configurations simulated with 
the Llama code, five of them being computed with two 


1 5 NR +(t) is monotonically increasing with t except for the 
datasets SXS:BBH:0178 and SXS:BBH:0065 where it is mono¬ 
tonically decreasing. 




FIG. 1. EOB/NR £ = m = 2 waveform comparison: q = 1, 
Xi = X 2 = 0 and ag(0.25) « —57.688 from Eq. |5|. The 
calibration was done using the 14-orbit simulation, following 
the same procedure as in Ref. [29] (that was also using the 
same dataset). The resulting EOB waveform is then checked 
for consistency against waveforms SXS:BBH:0066 (~ 28 or¬ 
bits, middle panel) and SXS:BBH:0002 (~ 32 orbits, bottom 
panel), that is the longest NR waveform currently available. 
The EOB/NR modulus and phase agreement is excellent all 
through inspiral, plunge, merger and ringdown. 


different resolutions. All but one of them have mass ratio 
q = 1. The dimensionless spin values that we consider 
are x = (±0.8, ±0.6, ±0.4, ±0.2, 0). The precise config¬ 
urations (and resolutions used) are listed in Table [II] 
The two nonspinning configurations were first presented 
in Ref. [261 : the remaining, equal-mass, spinning con- 
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figurations were first presented in Ref. [501 . They were 
here simulated at higher resolutions and with improved 
grid setups to enhance accuracy. For all Llama data the 
waveforms at future null infinity are estimated by using 
Cauchy-characteristic extraction (CCE) [3H 021 ■ 

The Llama simulations are relatively short (between 4 
and 8 orbits). We use them mainly for reaching one of 
the main aims of this paper, namely to check to what 
extent the EOB prediction for the energetics of the sys¬ 
tem (as measured by the gauge-invariant relation Eb(j) 
between the dimensionless binding energy and the angu¬ 
lar momentum) is consistent with the corresponding NR 
quantity. 

The longer SXS waveforms (between 12 and 32 orbits) 
not only were used to check the phasing performance 
of the EOB model, but also its energetics (modulo sub¬ 
tleties related to junk radiation discussed below, which 
are essentially absent when dealing with Llama data). 
The joined use of Llama and SXS waveform data al¬ 
lowed us to: (i) reach a reliable NR-calibration of the 
spinning EOB model of Paper I; (ii) perform interest¬ 
ing cross checks on energy and angular momentum at 
merger between NR data obtained with completely dif¬ 
ferent codes; and (iii) to compute the Eb(j) energetics 
for a large sample of spin-aligned SXS data (thereby im¬ 
proving on Ref. El where only two £/(/) curves were 
computed). 


III. IMPROVED CALIBRATION OF THE 
NONSPINNING EOBNR MODEL 

As mentioned in the Introduction, the spinning EOB 
model considered in this paper is essentially the one in¬ 
troduced in Paper I (to which we refer for our notation), 
in particular the waveform the waveform up to merger 
is defined by Eqs. (75)-(96) there. The only change in 
the theoretical frameworm concerns the definition of the 
EOB radial potential A. In the present work, we use of 
the full 4PN-accurate analytical knowledge of the EOB 
radial interaction (see below). In addition we extended 
(without any additional theoretical modification) the ap¬ 
plication of the new ringdown description of Ref. [14] 
to the full SXS datasets of Table [i] For the equal-mass, 
equal-spin configurations, the ringdown description is ex¬ 
actly as described in Ref. [14], including the use of the 
fitting coefficients listed in Table II there. On the other 
hand, for the remaining configurations in our Table [I] we 
do not make use of the latter fit, but rather we apply the 
ringdown modelization methodology of [14] separately to 
each dataset. We postpone to future work the construc¬ 
tion of a global analytical extrapolation of the latter, dis¬ 
crete, ringdown data. 

We use here the 5PNlog-accurate post-Newtonian ex¬ 
pansion of the orbital A™ function 

A™ (it c ) = 1 — 2u c + 2 uu\ + 

+ v[a%(v) + a^ e lnu c }ul + v[a c 6 {y) + ag° s lnuju®, (1) 


where 

= _M_ = mim 2 = g , . 

M (mi + m 2) 2 (1 + q) 2 

is the main deformation parameter of EOB theory, which 
varies between v = 0, in the large-mass-ratio limit 
(q = mi/m? 1) and v = 1/4 in the equal-mass case. 
The dimensionless gravitational potential u c is defined as 
u c = Mjr c in terms of the EOB centrifugal radius 


+ 5a 2 (r ) . (3) 

Here the (next-to-leading-order) correction of the Kerr 
parameter, Sa 2 (r ) [T2], is defined in Eq. (59) of Paper I. 
Contrary to Paper I, where we had phenomenologically 
fixed the 4PN coefficient a%{v) to the (^-independent) 
fiducial value a% = 23.5, we use here the exact, in¬ 
dependent, analytical expression of a^(v) obtained in 
Ref. [33] (see also [31]). We recall that the orbital 
EOB radial potential A or b is defined by Pade resumming 
A orb M as 


2 Ma 2 


Aorb(w c ; V, a c e) = PsIAorbK)]- (4) 

In view of the change in the analytical expression of A or b, 
our first task will be to provide a new calibration of the 
single, effective 5PN functional parameter a§(m) entering 
A or b- We perform this calibration by means of a sample 
of nonspinning waveforms. [This nonspinning-calibrated 
orbital potential will then be used as is in our spinning 
EOB model]. We use as calibrating waveforms eight SXS 
simulations with mass ratio q = (1,1.5, 2,3,4, 5, 6, 8) [the 
q = 9.989 one is used just as a cross check]. For q = 1 
we use the same ~ 14 orbits waveform that was used in 
the calibration procedure of Ref. [29] , We tune ae(n) so 
that the EOB and NR phasing agree (after a suitable 
alignment) within the NR phasing error at NR merger. 
Following the footsteps of Ref. [22] (and in particular 
the cross-check of the time-domain analysis with the Q u 
analysis) we got 

a%{v) = 3097.3 v 2 - 1330.6 v + 81.38. (5) 

Note that a%{v) varies between —57.69 for v = 0.25 
(equal-mass case) and —7.432 for q = 9.989 {v ~ 0.0827). 
Moreover the parameter AfNQc, giving the time-lag be¬ 
tween the peak of the pure orbital frequency fl 0 rb [see 
Eq. (100) of Paper I] and the NR/EOB matching point on 
the time axis (which is used to determine next-to-quasi- 
circular (NQC) corrections to the waveform) is fixed (as 
discussed in Sec. V of Paper I) as AtNQc(x = 0) = 1 M 
in the nonspinning caseQ Finally, as mentioned in the in¬ 
troduction, we use the new modelization of the ringdown 


2 By contrast, in the case of high, positive spins, \ — 0.85, it is 
useful to allow At^QC to depend on the spin (see below). 
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introduced in Ref. M- In this paper, we use the actual 
fitting coefficients as obtained applying the procedure of 
Ref. [H] to each single waveform. The performance of 
suitable “fits of fits” extending the validity of the proce¬ 
dure of Ref. [T3] outside the domain of known waveforms 
will be discussed elsewhere. 

Figure |T] (which plots the Regge-Wheeler-Zerilli nor¬ 
malized waveform ’F 22 = Rh^il V24) illustrates the 
EOB/NR phase agreement, for q = 1 when using: (i) 
the 14-orbit (up to merger) calibration waveform (top- 
panel); (ii) the ~ 28-orbit SXS:BBH:0066 waveform of 
the catalog; and (iii) the longest NR waveform available 
in the catalog, SXS:BBH:0002, that corresponds to about 
32 orbits up to merger. We stress that the actual cali¬ 
bration of a§(^) made use neither of SXS:BBH:0066 nor 
of SXS:BBH:0002, which are used here just to provide 
an independent check of the calibration procedure. The 
vertical lines in the plots highlight the alignment time in¬ 
tervals, corresponding to dimensionless GW frequencies 
(Mwi,Mwj{). To check consistency, we fixed Mwjj = 
0.045 for each data set. Then, we use Mujl = 0.035 
for the 14-orbit run, Mujl = 0.023 for SXS:BBH:0066 
and Muj l = 0.025 for SXS:BBH:0002. In the latter 
case, the phase difference between EOB and NR, accu¬ 
mulated up to NR merger, A<^‘( > ® NR = (() R ° g B — = 

0.073 rad is comparable to the NR uncertainty at merger, 
~ 0.066 rad (see Table]!]). Similarly, when using a%(v) as 
given by Eq. one finds that the EOB-NR phase dif¬ 
ference is comparable to (though in some cases slightly 
larger than) the NR uncertainty at merger for all 9 non¬ 
spinning data (compare columns seven and eight in Ta- 
ble|). 


IV. ENERGETICS OF NONSPINNING 
COALESCENCES: THE CHOICE T r , = 0 


A. Energetics with (nonspinning) Llama data 


Having used SXS phasing data to determine the rela¬ 
tion Eq. ([ 5 ]), let us now turn to discussing the energet¬ 
ics of the model, motivating, in particular, the choice 
of a vanishing radial component of radiation reaction 
that we made here. The choice T rt = 0, was intro¬ 
duced in Ref. [35] and used in several previous EOB 
works [261 [Ml EZ]) as well as in Paper I. As in previous 
works El 133 [M], the analysis of the energetics is done 
via the gauge-invariant relation between the dimension¬ 
less binding energy Eb and the dimensionless total angu¬ 
lar momentum, j, Eb(j). These quantities are computed 
as 


E b = 


H^ADM A£ r ad M 


■ _ ^ADM At/rad 

J ~ Mu 


( 6 ) 

(7) 


Here (A/^ DM , J^ DM ) denote the total, initial Arnowitt- 
Deser-Misner (ADM) mass-energy and angular momen- 



FIG. 2. Energetics for q = 1, nonspinning binary: com¬ 
parison between NR Llama data, the 4PN, Taylor-expanded, 
curve and two EOB curves, one with and the other without 
a radial part of the radiation reaction T r , ■ The (effective) 
choice J- r , = 0 displays the smallest (~ 10~ 4 ) discrepancy 
with NR data up to merger (indicated by colored markers). 
The inset focuses on the initial transient driven by junk radi¬ 
ation. 


turn of the system (including the contribution of the indi¬ 
vidual spins), (A£ ra( j, A^ad) denote the energy and an¬ 
gular momentum radiated in GWs, while M = mi + m 2 
and /i = TO 1 TO 2 /M, where mi and m 2 are the NR mea¬ 
sured initial Cliristodoulou masses. We recall that, while 
A£ ra d is obtained by time-integrating the square of the 
news (h), Aj7r a d is obtained by time-integrating a bilin¬ 
ear expression in the news and the strain h. To obtain the 
strain we start from the Fourier transform of ip 4 and use 
the method of [35] • The E h (j ) relation can be straightfor¬ 
wardly obtained from the results of the BBH simulations 
obtained with the Llama code [40] [44] , as we did in a pre¬ 
vious study |26j limited to nonspinning BBHs. On the 
other hand, when using SXS data, several subtleties have 
to be properly addressed in order to accurately estimate 
the Eb(j) relation (see Appendix for details). 

Figure [ 2 ] displays the Eb{j) relation in the q = 1, non¬ 
spinning case. It shows the triple comparison between 
the NR curve and: (i) the PN-expanded Eb(j) relation; 
(ii) the EOB model with a% given by Eq. (]5j) and with 
the choice T rt = 0; (iii)exactly the same EOB model (in 
particular, with the value of <Zg given by Eq. <§) but 
with the (nonresummed), 2.5PN accurate, expression of 
of Ref. [323 ■ Note that the latter expression was al¬ 
ready used in previous EOB work [22], though its effect 
on Eb(j) is investigated here for the first time. 

The general structure of the PN-expanded 
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FIG. 3. Energetics comparison, as in Fig. [2] but for q = 2, 
nonspinning, Llama NR data. The T r , = 0 curve displays the 
smallest discrepancy with NR data up to merger (indicated 
by colored markers). 


Eb(J',x I1X2) relation is (in the spinning case) 

= t ^. ( 8 ) 

71=2 

where £ = j — ( —-yq d- ~Xn ) denotes the dimension- 

\ m 2 to 1 ' J 

less orbital angular momentum, £ = L/([xM) and the 
coefficients c n (v, \) depend at most quadratically on the 
dimensionless spins %i = Si/ml, — S 2 /m The accu¬ 
racy we use in E™ is: 4PN in the nonspinning sector 03]; 
NNLO for the spin-orbit and NLO for the spin-spin part 
(see e.g. [ 00 ]). 

As shown in Fig . [ 2 ] (see especially the bottom panel, 
which displays A£f* NR (j) = E^ OB (j) - Ej/ R {j)) the 
use of a vanishing radial radiation reaction leads, from an 
effective point of view, to a better agreement between the 
EOB energetics and the NR one up to mergeij^] The NR 
and EOB merger points are indicated, respectively, by a 
filled circle (NR), by a filled square (EOB, = 0) and 
by an empty square (EOB, 7 ^ 0). By contrast, the 
PN curve shows the largest deviation from NR results, 
especially at low j' s. Finally, Fig. [3]shows that the same 
conclusions hold true also for the q = 2 Llama data. 

As discussed in Refs. S3 , there exists a coordinate 
gauge where F r , = 0. Further work j42j has shown that 


3 We found that a different calibration of a): is unable to displace 
the T rt 78 0 curve on top of the NR one 




FIG. 4. Comparison between EOB and SXS energetics for 
the q = 8 case. The shaded area indicates numerical uncer¬ 
tainty measured taking the difference between the highest and 
second highest resolution data. Note that the EOB merger 
state (red, filled, square marker) and the NR merger state 
(black, filled, circle) are barely distinguishable in the upper 
panel. 


the requirement of the vanishing of the Schott contribu¬ 
tion J Schott to the angular momentum required a specific, 
nonzero, value of T rt (which was then used in Ref. [25j 
because the vanishing of J Schott seemed apriori required 
by the definition of the azimuthal component T v of the 
radiation reaction within the EOB formalism). However, 
because of the various approximations made in defining 
T v in the EOB formalism, it is not actually required to 
choose the that is implied by a vanishing j Schott . 
Our present work experimentally shows that the condi¬ 
tion J- Tm = 0 is a simple and effective way to accomplish 
a (rather surprisingly) good agreement between the nu¬ 
merical and EOB (j ) curves. We leave to future work a 
theoretical study of why such a condition, in conjunction 
with the current approximate value of J rROB , happens to 
lead to rather small values of both J Schott and E^ chott 
(as shown in Figs. [2|3] above). 

B. Energetics with (nonspinning) SXS data 

We found that a similar excellent agreement among 
energetics holds, up to merger for all nonspinning SXS 
configurations at our disposal, which cover the mass-ratio 
range 1 < q < 9.989. [For a previous comparison of 
nonspinning merger characteristics see Ref. [05] ]. 

So far, the only published computation of the ener- 
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TABLE I. EOB/NR phasing comparison. The columns report: the number of the dataset; the name of the configuration in the 
SXS catalog; the mass ratio q = mi/m 2 ; the symmetric mass ratio v; the dimensionless spins xi and X+ the phase difference 
A0 eobnr = tp EOB — 0 NR computed at NR merger; the NR phase uncertainty at NR merger S()>^ g (when available) measured 
taking the difference between the two highest resolution levels (see text); the maximum value of the unfaithfulness F = 1 — F 
as per Eq. ([22]) . The A</> BOBNR ’s in brackets for \i = X 2 > +0.85 were obtained using Eq. plj) for Af NQC (x)- 


# 

Name 

N orbits 

q 

V 

Xi 

X2 

a 1EOBNR r„_ ji 

A <p mrg [rad] 

<5<(>m?g [rad] 

max(E) 

1 

SXS:BBH:none 

14 

1 

0.25 

0.0 

0.0 

-0.016 


0.00087 

2 

SXS:BBH:0066 

28 

1 

0.25 

0.0 

0.0 

+0.010 


0.00068 

3 

SXS:BBH:0002 

32.42 

1 

0.25 

0.0 

0.0 

+0.073 

0.066 

0.00101 

4 

SXS:BBH:0007 

29.09 

1.5 

0.24 

0 

0 

+0.05 

0.018 

0.00201 

5 

SXS:BBH:0169 

15.68 

2 

0.2 

0 

0 

-0.15 

0.02 

0.00045 

6 

SXS:BBH:0030 

18.22 

3 

0.1875 

0 

0 

-0.074 

0.087 

0.00035 

7 

SXS:BBH:0167 

15.59 

4 

0.16 

0 

0 

-0.059 

0.52 

0.00035 

8 

SXS:BBH:0056 

28.81 

5 

0.138 

0 

0 

-0.089 

0.44 

0.00038 

9 

SXS:BBH:0166 

21.56 

6 

0.1224 

0 

0 

-0.198 


0.00037 

10 

SXS:BBH:0063 

25.83 

8 

0.0987 

0 

0 

-0.453 

1.01 

0.00292 

11 

SXS:BBH:0185 

24.91 

9.98911 

0.0827 

0 

0 

-0.0051 

0.376 

0.00066 

12 

SXS:BBH:0004 

30.19 

1 

0.25 

-0.50 

0.0 

-0.017 

0.068 

0.00403 

13 

SXS:BBH:0005 

30.19 

1 

0.25 

+0.50 

0.0 

+0.08 

0.28 

0.00052 

14 

SXS:BBH:0156 

12.42 

1 

0.25 

-0.95 

-0.95 

+0.32 

2.17 

0.00058 

15 

SXS:BBH:0159 

12.67 

1 

0.25 

-0.90 

-0.90 

+0.06 

0.38 

0.00047 

16 

SXS:BBH:0154 

13.24 

1 

0.25 

-0.80 

-0.80 

+0.11 


0.00044 

17 

SXS:BBH:0151 

14.48 

1 

0.25 

-0.60 

-0.60 

-0.049 

0.14 

0.00042 

18 

SXS:BBH:0148 

15.49 

1 

0.25 

-0.44 

-0.44 

+0.14 

0.72 

0.00043 

19 

SXS:BBH:0149 

17.12 

1 

0.25 

-0.20 

-0.20 

+0.45 

0.90 

0.00085 

20 

SXS:BBH:0150 

19.82 

1 

0.25 

+0.20 

+0.20 

+0.94 

0.99 

0.00275 

21 

SXS:BBH:0152 

22.64 

1 

0.25 

+0.60 

+0.60 

+0.01 

0.36 

0.00068 

22 

SXS:BBH:0155 

24.09 

1 

0.25 

+0.80 

+0.80 

-0.39 

0.26 

0.00110 

23 

SXS:BBH:0153 

24.49 

1 

0.25 

+0.85 

+0.85 

+0.06 


0.00059 

24 

SXS:BBH:0160 

24.83 

1 

0.25 

+0.90 

+0.90 

+0.41 (+0.41) 

0.80 

0.00117 

25 

SXS:BBH:0157 

25.15 

1 

0.25 

+0.95 

+0.95 

+0.37 (+0.83) 

1.18 

0.00295 

26 

SXS:BBH:0158 

25.27 

1 

0.25 

+0.97 

+0.97 

+0.37 (+0.49) 

1.26 

0.00325 

27 

SXS:BBH:0172 

25.35 

1 

0.25 

+0.98 

+0.98 

+0.99 (+0.46) 

2.02 

0.00422 

28 

SXS:BBH:0177 

25.40 

1 

0.25 

+0.99 

+0.99 

+0.22 (+0.48) 

0.40 

0.00507 

29 

SXS:BBH:0178 

25.43 

1 

0.25 

+0.994 

+0.994 

+0.24 (+0.23) 

-0.53 

0.00506 

30 

SXS:BBH:0013 

23.75 

1.5 

0.24 

+0.5 

0 

+0.31 


0.00058 

31 

SXS:BBH:0014 

22.63 

1.5 

0.24 

-0.5 

0 

-0.15 

0.15 

0.00046 

32 

SXS:BBH:0162 

18.61 

2 

0.2 

+0.6 

0 

-0.20 

0.71 

0.00027 

33 

SXS:BBH:0036 

31.72 

3 

0.1875 

-0.5 

0 

+0.08 

0.065 

0.00040 

34 

SXS:BBH:0031 

21.89 

3 

0.1875 

+0.5 

0 

+0.12 

0.034 

0.00023 

35 

SXS:BBH:0047 

22.72 

3 

0.1875 

+0.5 

+0.5 

-0.034 


0.00030 

36 

SXS:BBH:0046 

14.39 

3 

0.1875 

-0.5 

-0.5 

+0.36 


0.00054 

37 

SXS:BBH:0110 

24.24 

5 

0.138 

+0.5 

0 

+0.24 


0.00016 

38 

SXS:BBH:0060 

23.17 

5 

0.138 

-0.5 

0 

+0.21 

0.8 

0.00034 

39 

SXS:BBH:0064 

19.16 

8 

0.0987 

-0.5 

0 

+0.026 

0.8 

0.00042 

40 

SXS:BBH:0065 

33.97 

8 

0.0987 

+0.5 

0 

+1.33 

-3.0 

0.00040 


getics from SXS data was presented in Ref. m- it 
was limited to two, quasi extremal, spinning configura¬ 
tions with spin either aligned or antialigned with the or¬ 
bital angular momentum, with Xi = Xi = ^0.95 and 


Xi = X 2 = +0.98. 

We present here the first systematic computation of 
E b (j) curves from SXS data covering a large sample 
of nonprecessing, spinning (and nonspinning) configura- 
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TABLE II. BBH configurations, evolved using the Llama code, used to compute the gauge-invariant relation between energy 
and angular momentum. From left to right: name of the run, number of orbits up to merger, mass ratio q = mi/m 2 , symmetric 
mass ratio v, dimensionless spins of the two black holes, (xi> X 2 ); the initial ADM mass M^ dm ; initial total angular momentum 
Jadm- 


Name 

N orbits 

Q 

V 

Xi 

X2 

A?adm 

✓7-0 

‘-'ADM 

ql_sO 

8.09 

1 

0.25 

0.0 

0.0 

0.99051968 

0.99325600 

q2_s0 

6.70 

2 

0.2 

0.0 

0.0 

0.990898 

0.85599600 

ql_s-8D10_h96 

3.83 

1 

0.25 

-0.8 

-0.8 

0.989412 

0.61736 

ql_s-8D10_hll52 

3.84 

1 

0.25 

-0.8 

-0.8 

0.989412 

0.61736 

ql_s-6D12_h64 

8.25 

1 

0.25 

-0.6 

-0.6 

0.986161 

0.7552392 

ql_s-6D12_h512 

8.23 

1 

0.25 

-0.6 

-0.6 

0.986161 

0.7552392 

ql_s-4DKLh64 

5.67 

1 

0.25 

-0.4 

-0.4 

0.990138 

0.791588 

ql_s-2D10_h64 

5.88 

1 

0.25 

-0.2 

-0.2 

0.989941 

0.877499 

ql_s-2D10_h512 

5.88 

1 

0.25 

-0.2 

-0.2 

0.989941 

0.877499 

ql_s2D8_h64 

4.93 

1 

0.25 

+0.2 

+0.2 

0.987658 

0.984992 

ql_s2D8_h512 

4.93 

1 

0.25 

+0.2 

+0.2 

0.987658 

0.984992 

ql_s4D8_h64 

5.67 

1 

0.25 

+0.4 

+0.4 

0.987619 

1.072176 

ql_s6D8_h64 

6.38 

1 

0.25 

+0.6 

+0.6 

0.987649 

1.160296 

ql_s6D8_h512 

6.48 

1 

0.25 

+0.6 

+0.6 

0.987649 

1.160296 

ql_s8D8_h64 

7.24 

1 

0.25 

+0.8 

+0.8 

0.987744 

1.249296 


x nr 4 



FIG. 5. EOB-NR differences between Eb(j) curves for all but 
one (q = 1.5) nonspinning SXS data of Table]!] The difference 
is compatible with NR uncertainties up to merger (~ 10“ 4 , 
see Fig. |4] for the illustrative q = 8 case). 


tions. We computed Ef,(j) curves from the nonspinning 
binaries in Table [I] applying the same standard procedure 
we used for Llama data. Note, however, that, while the 
time-integration of energy and angular momentum losses 
from Llama data could be meaningfully performed start¬ 
ing at the beginning of the simulation, in the case of SXS 
data, we found necessary to correct the time-integration 
of the early “junk-radiation” losses by adding a vecto- 



FIG. 6. Binding energy and angular momentum at merger 
(i.e. at the maximum of |ft. 22 1) for nonspinning binaries versus 
the symmetric mass ratio v. The plot shows: (i) the highly 
accurate numerical agreement between SXS and Llama final 
configurations for q — 1 and q = 2; (ii) the good numerical 
consistency between EOB predictions and actual NR states. 


rial shift (A j°, AEjj) in the (J,Ei,) plane. This shift was 
essentially determined so as to minimize the EOB/NR 
difference AE^ OBNR (j) during the early inspiral. All 
technical details of our computation are explained in the 
Appendix [see in particular Table [TTT] that lists the values 
of the vectorial shifts (A j°, AEjj) we used]. 














The main results of our SXS-data analysis are: (i) for 
q = 1 and q = 2, the energetics computed from SXS data 
confirms the results based on Llama data given above in 
Figs. H (ii) for the other values of q, 3 < q < 9.989, not 
covered by Llama data, we again find an excellent con¬ 
sistency between the energetics of our EOB model and 
the SXS NR one; (iii) it is in particular remarkable that 
the binding energy and angular momentum at merger 
(E™ rg , j“ rg ) provided by our EOB model are in accurate 
agreement with the NR one R These results are exem¬ 
plified in Figs. 00 and Fig. [6J 

For illustrative purposes, Fig. [4] explicitly shows the 
q = 8 case, dataset SXS:BBH:0063. The bottom panel 
shows the EOB-NR difference. The vectorial shift in the 
(J, E b ) plane used here is A j° = —4 x 10~ 3 , A E b = 
-7.7709 x 10 -4 (so that E b (j) = E T b aw (j - A j°) + A E$. 
As discussed in detail in the Appendix these (quite small) 
shifts are determined through the following two-step pro¬ 
cedure. Step 1: A j° is determined first, by requiring 
that the difference between the EOB and NR dimen¬ 
sionless orbital frequencies (defined as 11 = dE b /dj) 


AfiEOBNR (i) = QEOB (j) _^NR (j) osci n ates around zero 
on the largest possible j-interval. With a so chosen Aj 0 , 
one finds that Afi EOB ^ R oscillates between ±5 x 10 -5 
over the interval 3.6 < j < 4.2. Step 2: the constant 
AE b is then determined by minimizing the difference 
AFfOBNR(j) == E BOB {j) - El™{j - A j°), which is es¬ 
sentially constant on a large j interval, during the early 
inspiral. After adding these small shifts to the raw data, 
one finds that AifE 0 ® 1411 ^) is of the order of ±10 -5 for 


3.5 < j < 4.2. 

The shaded area in the bottom panel of Fig. [4] indi¬ 
cates the numerical uncertainty, as estimated by taking 
the difference between the (raw) E b (j) curves obtained 
from the highest (Lev5) and second-highest (Lev4) reso¬ 
lution data present in the SXS catalog. In order to get a 
conservative uncertainty estimate, we do not apply any 
vectorial shift to the raw data. Let us emphasize the re¬ 
markable agreement between EOB and NR energetics at 
merger. The difference at the NR merger point is of order 
4 x 10 —4 , which is barely visible in the top panel of Fig.0 
(and it is approximately twice the estimated uncertainty 
there). 

For the other nonspinning configurations in Table |T| 
one obtains curves quite analogous to the q = 8 case. The 
relevant information is displayed in Fig. [5] that collects 
the A£’ eobnr (j) = Ef OB (j) — E^ R (j) differences for all 
nonspinning data but one (the q = 1.5 one). Note that 
the end points indicated by square markers in Fig. [5j 
record the values of the function AE eobnr (j) at the 
corresponding EOB merger values of j, j E ° B . On the 
other hand, the stars on the axaxis of the figure mark the 
corresponding NR merger values of j, _)^ R g . 


For q = 1, we use here the SXS:BBH:0066 dataset with 
~ 28 orbits. The vectorial shifts (Aj b , AE b ) we use are 
listed in Table m in the Appendix. For all configura¬ 
tions, one finds that the EOB/NR agreement is of the 
order 10~ 5 at the beginning of the inspiral and grows 
only up to a few parts in 10 ~ 4 around merger. These dif¬ 
ferences are comparable to (though slightly larger than) 
the corresponding error bars computed, as above, taking 
the difference between the two highest resolutions avail¬ 
able. 

In Fig. [ 6 ] we complement the above results by focus¬ 
ing on the EOB/NR comparison for merger characteris¬ 
tics, (jrnrgj-E/nrg) as a function of v. The figure shows 
that the EOB/NR disagreement increases with v. from 
a difference of 3.07 x ICC 4 for E“ rg and 1.7 x ICC 4 for 
j mr g for q = 9.989 [y = 0.0827) the EOB model ends 
up slightly overestimating the NR values in the equal- 
mass case, v —> 0.25. In particular, for v = 0.25 we 
see that the EOB prediction for E™ 6 is larger than 
the NR one by 4.9 x ICC 3 and that the EOB predic¬ 
tion for j mrg is larger than the NR one by 3.1 x ICC 2 . 
[These differences are larger than those apparent in the 
corresponding endpoints of Fig. [5] because we are now 
computing E BOB {j B< ? B ) - E™(j" B ) while the vertical 
displacement of the endpoints of Fig. [5] corresponded to 
^h EOB (jm? g B ) ^ ^h NR 0 m? g B )]- In the same figure we also 
display Llama merger data, denoted by asterisks. We 
find that Llama and SXS data are fully consistent within 
the ICC 4 level. Finally, the ^-dependence of the SXS 
NR E™ rg can be fitted by means of a second order poly¬ 
nomial, while j mrg (i/) is well fitted by a straight line, 
namely 

E ™| s = -0.062951 - 0.176949i/ - 0.216264i/ 2 , (9) 

j™xs = 3.406583 - 2.421052^. (10) 

V. NR CALIBRATION OF THE EOB 
SPIN-ORBIT COUPLING 

Let us now turn to discussing the numerical comple¬ 
tion of the EOB spinning model. As mentioned above, 
the spinning sector of the model we use here is exactly 
the same as described in Paper I. In the latter work it 
was shown that it was possible to obtain a good EOBNR 
phasing agreement by tuning a single effective functional 
parameter C 3 (x) that was entering the spin-orbit coupling 
at next-to-next-to-next leading ordeijj The best values 
of the parameter C 3 , for various spin values, are listed 
in Table I of Paper I. Since the orbital A m b function has 
changed, in the present work, because of (i) the use of the 
complete analytical 4PN information and (ii) the conse¬ 
quently modified functional dependence of a®(^) given 


4 As mentioned above, we recall that the instant of the merger 

is defined, both in EOB and NR, as the time when the of the 5 The analysis of Paper I was limited to the equal-mass, equal-spin 
modulus of the l = m = 2 waveform reaches its maximum. case. 
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by Eq. ([5]), we now need to look for a new determination 
of C 3 (x). In doing so, we also consider here more NR 
simulations than in Paper I, notably taking into account 
all available nonprecessing data in the SXS catalog. The 
spinning SXS configurations we use are listed in Table |T] 
the mass ratio varies in the range 1 < q < 8 and there are 
several configurations where only one of the two black- 
holes is spinning. A priori, one expects the NNNLO ef¬ 
fective parameter C 3 to be a function of both the mass 
ratio and the spins of the binary. The determination of C 3 
is done in two steps. In a first step, we separately consid¬ 
ered each binary configuration and determine a prelimi¬ 
nary best value of C 3 by minimizing the EOB-NR phase 
difference, after alignment in the inspiral phase, so to be 
compatible with (and typically smaller than) the NR un¬ 
certainty. This procedure is rather straightforward, as it 
is just a one parameter search. In a second step, we look 
for a global, analytical representation that approximately 
reproduces the latter preliminary best values of C 3 as a 
function of symmetric mass ratio and spins. We found 
that one can represent, with sufficient accuracy, the val¬ 
ues (determined by minimizing the EOB/NR phase dif¬ 
ference as explained above) of C 3 for the entire sample of 
configurations listed in Table [TJ by means of the following 
simple functional relation 


c 3 (ai,a 2 ,i / ) =Po 


1 + ni(di + a 2 ) + n 2 (ai + a 2 ) 2 
1 + d\ (a 2 + a 2 ) 


+ (piv + p 2 v 2 + P 2 V 3 )(ai + a 2 )V 1 - 4 v 


+ Pi{ai - a 2 )v 2 , (11) 

where 

po = 44.786477, (12) 

ni = -1.879350, (13) 

n 2 = 0.894242, (14) 

di = -0.797702, (15) 

pi = 1222.36, (16) 

p 2 = -12764.4, (17) 

p 3 = 36689.6, (18) 

p 4 = -358.086. (19) 


and where we found convenient to introduce the spin 
quantities di j2 = A' 1 j 2 Xi, 2 , with X\ 2 = toi j 2 /M, and 
M = m\ + m 2 . With our convention that mi > m 2 , in 
terms of the symmetric mass ratio v we have 

Ad = \ (1 + Vl^v) , (20) 

and X 2 = 1 — Ai. The terms in Eq. © that vanish in 
the equal-mass, equal-spin case were chosen, for simplic¬ 
ity, to be linear in the spins. Similarly, the polynomial 
dependence in v was found necessary to properly fit the 
values of C 3 yielding a good NR/EOB phasing agreement 
for q = 8 , (xi,X 2 ) = (+0.5,0), SXS:BBH:0065 configu¬ 
ration. In the equal-mass, equal-spin (aligned, or anti¬ 
aligned with the orbital angular momentum) case, one 



FIG. 7. Quasi-linear behavior of the NNNLO effective 
spin-orbit coefficient 03 ( 01 , 02 , v), Eq. (Ill, in the equal-mass 
( 1 / = 0.25), equal-spin ai = 02 = X case. The tuning 
of this single dynamical parameter allows to get an excel¬ 
lent EOB/NR phasing agreement throughout inspiral, plunge, 
merger and ringdown. 


finds that the dependence of C 3 on the spin is nearly lin¬ 
ear (see Fig. [ 7 ]). 

Finally, let us mention that, as already discussed in 
Paper I (see Table I there) we found it necessary to 
flex the simple choice AtNQC = 1M, uniformly used 
in the nonspinning case, so as to allow it to depend 
on spin for large, positive, spins. Actually, the only 
available simulations where we found the need of flexing 
A^qc a re the six equal-mass, equal-spin configurations 
with x = Xi = X 2 > Xo = +0.85. In practice, for 
spins x = {+0.90, +0.95, +0.97, +0.98, +0.99, +0.994} 
we found as good choices AfNQC = 
{0.2,—1.2,—1.7,—2.0,—3.0,—3.2}, respectively. In 
the EOB numerical evolution we use a time-resolution 
At EOB = 0.1 M, and Af N QC is chosen as an integer 
multiple of At EOB . The values of AfNQC listed above 
are accurately fitted by using a (1,1) Pade approximant: 


AtNQC (x) 


1 + 7+(x - Xo) 
i + <Mx-xo) 


( 21 ) 


with ni = —16.06288 and d\ = —4.04266 and xo = 0.85. 

The quality of the fit yielded by Eq. © (together 
with the discrete values of AtNQC listed above) is quan¬ 
titatively assessed by measuring the EOB-NR phase dif¬ 
ference at NR merger after having aligned (in time and 
phase) the EOB waveforms to the NR waveform during 
the early inspiral. Such differences are listed as A</) E ° BNR 
in Table[I] The same table also clearly illustrates the com¬ 
patibility of the EOBNR model with the numerical phase 
uncertainties at merger all over the waveform sam¬ 

ple considered. The use, in addition to Eq. (11), of the fit 
(21) slightly worsens A</) E ° BNR as indicated by the val- 
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FIG. 9. Unfaithfulness of the t = m = 2 EOB waveforms with 
respect to the NR ones for all BBH datesets of Table [I] The 
largest values of F corresponds to the quasi-extremal spinning 
configurations SXS:BBH:0177 and SXS:BBH:0178 in Table [ij 
The labels of the configurations use the notation (q,x 11 X 2 ); 
they are highlighted in color so as to ease a direct comparison 
with Fig. 1 of Ref. El- 


model presented here, and to give a clearer physical 
meaning to the phase differences quoted above, we also 
measured the agreement between the EOB waveforms 
and all the available NR ones by computing the EOB/NR 
unfaithfulness (as a function of the total mass M) 


FIG. 8. Illustrative EOBNR time-domain phas¬ 

ing comparison for quasi-extremally spinning binaries: 
( 9 , Xi, ») = (!,+0.98,+0.98), dataset SXS:BBH:0172 
(top panel), and (q, yi, x' 2 ) = (1, +0.994, +0.994), dataset 
SXS:BBH:0178 (bottom panel). In both cases, the EOBNR 
difference at merger (dashed vertical line) is within the 
corresponding NR uncertainty, that is ~ 2 rad for y — +0.98 
and ~ 0.53 rad fory = +0.994 (see Table [j|. 


ues in parantheses in Table |T] Anyway, both values are 
approximately within (half of) the numerical error bar. 

Let us emphasize that the EOB/NR agreement re¬ 
mains excellent for the near-extremal spinning cases xi = 
X '2 = +0.98 and Xi = X 2 = +0.994. Figure [t] shows the 
EOB/NR comparison obtained by performing the usual 
time-domain comparison when aligning the waveforms 
on the frequency interval (Mwi, Muir) corresponding to 
the time-interval indicated by the two vertical dashed 
line in the plot. Analogous (or better) plots are found 
for all other configurations. As an indicator of this good 
EOB/NR agreement we just give in Table [I] the values 
of the phase difference at merger. 

To further demonstrate the high-quality of the EOB 


F(M ) = 1 - max ^ 2 ° B > 


( 22 ) 


where t 0 and </ 0 are the initial time and phase, ||h|| = 
y/ (h, h), and the inner product between two waveforms 
is defined as (/u, h 2 ) = 45ft J(nr (M) /S n {f) df , 

where S n {f) is the zero-detuned, high-power noise spec¬ 
tral density of advanced LIGO .47] and //J?(M) = 
is the starting frequency of the NR waveform (af¬ 
ter the junk radiation initial transient). Both EOB and 
NR waveforms are tapered in the time-domain so as to 
reduce high-frequency oscillations in the corresponding 
Fourier transforms. The procedure of tapering the wave¬ 
forms is the same followed by Ref. EL that introduced 
a different EOB model calibrated to NR, called SEOB- 
NRv2, and performed the same unfaithfulness analysis 
we are doing here. Figure [9] shows the so computed un¬ 
faithfulness as a function of the total mass of the binary 
for all configurations we considered. The maximum of 
F(M) is also listed, for convenience, in the last column 
of Table U One sees that for most of all considered con¬ 
figurations the value of F always stays one order of mag¬ 
nitude below the reference value of 1% (actually, most 
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FIG. 10. EOBNR time-domain phasing comparison for 
SXS:BBH:0065 configuration, q = 8, (XI 1 X 2 ) = (+0.5,0). 
The EOBNR difference at merger (dashed vertical line) is 
compatible with the corresponding NR uncertainty ~ —3 rad 
(see Table [I|. The eccentricity (~ 10~ 3 ) modulation in the 
phase difference is rather visible. 


configurations range between 0.1% and 0.01%) as the to¬ 
tal mass of the binary ranges from 20 to 200 Mq. Such 
a waveform quality implies a negligible loss in event rate 
due to the modeling uncertainty within the frequency 
range / > Pj(M). Note,however, that the NR wave¬ 
forms do not cover the entire frequency band of the de¬ 
tector when M < 1OOM 0 . For example, the longest 
NR waveform available, SXS:BBH:0002, has /^(M) ss 
900 (Mq/M) Hz, so that it starts at a frequency < 10 Hz 
when M > 90 Mg (and < 20 Hz when M > 45 M 0 ). The 
corresponding value of the unfaithfulness for M = 9OM 0 
is F = 8.54 x 10 -4 . For the highest spinning waveform, 
SXS:BBH:0178, we have /££(M) « 13OO(M 0 /M) Hz 
so that we need M > 13OM 0 to reach = 10 Hz. 
However, several recent works have shown that the EOB 
formalism provided the best available description of GW 
phasing even during early inspiral |48H50| . In absence 
of longer NR waveforms able to cover the entire fre¬ 
quency range of the detector for M > 20M 0 , the val¬ 
ues of F(M ) displayed in Fig. [ 9 ] do not have a direct 
data-analysis meaning over the full plotted mass range 
20 Mq < M < 200M 0 . We, however, expect, notably in 
view of Ref. [49], that the values of F will not significantly 
degrade when using longer NR waveforms. 

Figure [9] highlights in color the same particular con¬ 
figurations that were highlighted in Fig. 1 of Ref. EH], 
so as to prompt an easy and direct comparison. In ad¬ 
dition, Fig. [9] also shows, in purple, the configuration 
(+0.994, +0.994) that was not available when Ref. |lT] 
was written. Note that the worst global unfaithfulness, of 
the order of max F ss 0.005, corresponds to the two quasi¬ 
extremal spinning cases Xi = X 2 = (+0.99,+0.994), see 


Table [TJ A direct comparison with Fig. 1 of Ref. QTj al¬ 
lows us to make the following observations: (i) the config¬ 
uration (q,xI1X2) = (1,+0.6,+0.6), SXS:BBH:0152, de¬ 
livers, within our EOB model, F(M) < 10 -3 , all over the 
total mass range considered. The corresponding curve in 
Fig. 1 of Ref. mi was starting around F ~ 4.5 x 10 3 
for M = 20 Mq, then increasing up to F ~ 10 -2 for 
M ~ 50 Mq before decreasing again down to 2.5 x 10~ 3 
for M = 200 Mq. Similary, the other rather extreme case 
(q,XUX 2 ) = (8, +0.5, 0.6), SXS:BBH:0065, yields here 
an unfaithfulness of order 4 x 10~ 4 at M = 20 M 0 , that 
then decreases by a factor two for larger total mass. By 
contrast, it is interesting to note that the EOB/NR un¬ 
faithfulness computed for the same configuration using 
SEOBNRv2 increases with M from « 0.3% to « 1% (see 
Fig. 1 of Ref. HU). 

Note that the, numerical phase uncertainty on this 
dataset is as large as —3 rad at merger, so that c 3 
was (conservatively) calibrated so as to yield a 1.3 rad 
EOB/NR phase difference at merger (see Fig. 10). We 
advocate new simulations with reduced error bars for this 
configuration to firm up (and possibly improve) the cur¬ 
rent EOB calibration. 

Globally, the results collected in Fig. [9] show that our 
model quantitatively improves upon existing results. We 
shall discuss more these and other aspects of our unfaith¬ 
fulness comparison in the Conclusions (see in particular 
Fig. 21 and related discussion there). 


VI. ENERGETICS FOR SPINNING 
COALESCENCES 

A. Energetics of spinning Llama data 


Let us finally discuss the energetics of spinning coa¬ 
lescences yielded by our newly calibrated EOB model. 
We start doing this with Llama data and we will cross 
check our results with SXS data in the next Section. Fig¬ 
ure [ll contrasts the NR and EOB curves with 


Xi = X 2 = (±0.2, ±0.4, ±0.6), with the EOB-NR dif¬ 
ference AEf,(j) shown in each bottom subpanel. As be¬ 
fore, the EOB (red) and NR (black) mergers are indi¬ 
cated by markers. One sees that the differences are of 
the order of 10 -4 (or less) during the inspiral, to grow 
up to approximately the 10 -3 level around merger. One 
also notices that the disagreement between NR and EOB 
merger quantities depends on the configuration. In sev¬ 
eral representative cases we have indicated in Fig. El an 
estimate of the NR uncertainty on the energetics. The 
latter estimate was obtained by first taking the differ¬ 
ence between the Eb(j) curves computed from the highest 
(As = 0.512) and second highest (As = 0.64) resolution 
at our disposal and then Richardson extrapolating it as¬ 
suming some convergence order. The three colored areas 
around the difference in the bottom panels of Fig. El dis¬ 
play three different estimates of this NR uncertainties, 
as obtained by: (i) conservatively assuming 2nd order 
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FIG. 11. Energy versus angular momentum curves for spinning equal-mass binaries with spin x either aligned or anti-aligned 
to the orbital angular momentum, as obtained from Llama data. Good mutual consistency EOB/NR is found. The shaded 
area indicates an error bar estimate obtained by Richardson extrapolating three resolutions. 


convergence (lighter area); (ii) assuming 4th order con¬ 
vergence; and (iii) assuming 8th order convergence. The 
most conservative estimate (light grey), may eventually 
dominate the error budget in the limit of infinite reso¬ 
lution because of the presence of second-order finite dif¬ 
ference operators in the numerical infrastructure. This 
gives bounds of the order ±10” 3 around merger. Such a 
10~ 3 level is compatible with the EOB/NR differences at 
merger that we find for all other configurations. Other 
sources of uncertainty on the NR data such as: (i) the 
conversion from to h(. m (39J; (ii) the CCE error due 
to the choice of the finite-radius worltube where the CCE 
is started 123 E]; or (iii) the effect of higher multipoles, 
that contribute at the 10 -4 level, are subdominant with 
respect to the resolution uncertainty. 


1. High spins: ,\i = X 2 = ±0.8 (Llama data) 


The higher spin values x = Xi = Xi = ±0.8 deserve 
a separate discussion. First of all, inspecting the Eb(j) 


curves one sees that, contrary to the previous cases, the 
junk radiation transient is such that the “raw” NR curves 
(depicted as black, dashed lines in Fig. 12) stand visibly 
below the EOB prediction (for both signs of spin), see in¬ 
sets of the figure. This phenomenon is analogous (though 
smaller in magnitude) to what we found in nonspinning 
SpEC data. As mentioned above, it suggests that the NR 
Eb(j) curve should be corrected by an additional shift 
vector (A j°, A E®) so to have consistency with the EOB 
prediction at large values of j. The phenomenon found 
here with Llama data is probably due to the fact that 
the higher multipoles (l > 4) of the junk-radiation are 
not resolved with sufficient accuracy at the spatial reso¬ 
lution we can afford (see below). The error bars, com¬ 
puted as above (though only for one configuration) are 
also rather large and essentially constant over the entire 
j-range. Assuming the conservative second-order conver¬ 
gence order, one sees that the observed energy difference 
is mostly compatible with the error bar, remaining rather 
flat up to j ~ 2. We then apply the same technique dis¬ 
cussed above to determine (Aj°,AE 3 ), i.e., we inspect 
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FIG. 12. Same as Fig. [TT] (still for Llama data), but for the 
cases % — ±0.8. The black curves are obtained after shifting 
the raw data (dashed line) with the vector (Aj°,AE°) so as 
to recover a good EOB/NR consistency for large values of j. 
The Taylor-expanded PN curve is also added for comparison. 
The insets zoom on the initial (junk-related) drop in the E^(j) 
curve. 


the EOB/NR differences Afi EOBNR (j) and A.E EOBNR (j) 
and determine the vectorial shift so as to have them as 
flat and as small as possible in the first part of the inspi¬ 
ral (see Appendix for a fuller discussion). We choose here 
(Aj°, AE°) = (—1.3x 10~ 3 , 7.8x 10~ 4 ) for x = -0.8 and 
(A j°, A El) = (-1.3 x 10~ 2 , 7 x 10" 4 ) for X = ±0.8. Fig¬ 
ure 13 (that refers to the x = ±0.8 case) illustrates the 
effect that a good choice of Aj° has in centering around 




FIG. 13. Comparison between the analytical and numerical 
dimensionless orbital frequencies Cl for the (Llama) case xi = 
X 2 = +0.80 in Fig. |12| The bottom panel shows how the 
difference Af2 EOBNR ( i j) = 0 EOB (j) — 0 NR (j) can be made 
to oscillate around zero by a proper choice of the angular 
momentum shift A j°. 


a zero averaged value the (eccentricity driven) oscilla¬ 
tion in AO EOBNR (j). Figure 12 (bottom panels) shows 


how the use of the vectorial shifts substantially reduces 
the difference Ai? EOBNR , making it compatible with the 
4th-order extrapolation error bar. 

Let us finally comment on the possible origin of the 
excessive initial drop in the E/(j) curve (visible in the 
insets of Fig. 12) yielded by the junk-radiation transient. 


We think it has to do with the under-resolution of the 
junk-radiation high-multipoles (4 < £ < 8) for the fol¬ 
lowing reason: we noted that the energy loss during the 
junk-transient yielded by the modes with 5 < £ < 8 
is approximately as large as that for 2 < £ < 4 for 
Xi = X 2 = ±0.8; on the contrary, for the other con¬ 
figurations one always observes a sort of convergence of 
the junk radiation losses, with the higher multipoles con¬ 
tributing vertical drops that are progressively smaller 
than those due to the leading order modes. This makes 
us suspect that if it were possible to improve the ac¬ 
curacy of the higher-modes junk radiation one would be 
able to find the same straightforward (i.e., without shifts) 
EOB/NR consistency as for the case —0.6 < X < +0.6. 
A more refined analysis of these issues (either employ¬ 
ing higher-resolution and/or starting the system at larger 
separations, so as to reduce the magnitude of the junk) 
is left to future work. 

For the moment, to gain more confidence in our analyt¬ 
ical model, as well as to gauge the evidence of systematic 
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uncertainties in the NR data, we turn to explore the en¬ 
ergetics of spinning binaries computed using SXS data. 


B. Energetics of SXS spinning data 


1. Equal-mass, spinning binaries: —0.8 < x < +0.8 


One of the crucial outcomes of the Llama/EOB com¬ 
parison seen above is that, even when the separations 
are relatively small, the NR energetics, represented by 
Eb{j), quickly settles down to the EOB one after the ini¬ 
tial junk radiation transient. This was clearly the case for 
—0.6 < X < +0.6. For x = ±0.8, we recovered approxi¬ 
mately the same behavior after properly fixing a suitable 
vector shift (A j°, A E°). Such a “convergence” of the NR 
Eb(j) curve to the EOB one after the initial transition, 
would, a priori, be expected to emerge even more clearly 
from any of the SXS dataset at our disposal, since the 
initial separation between the black holes is much larger, 
so that the EOB curve should give an even better approx¬ 
imation to the NR one. However, in practice, as detailed 
in the Appendix, the SXS junk radiation is often rather 
large and tends to produce unphysical effects in the Eb(j) 
curve. In some cases, the “raw” Eb(j) NR curve is shifted 
in a region of the (j, Eb) plane which is clearly incompat¬ 
ible with the EOB prediction even in the early inspiral 
(see in particular Fig. 25 in the Appendix). As above, 
the problem is overcome by introducing a suitable shift 
vector (Aj°, A E°). More technical details, as well as the 
values of the shifts, are given in the Appendix. 

Figure [l4] illustrates the Eb(j) comparison between the 
EOB predictions and the SXS data. The figure dis¬ 
plays just the restricted sample of SXS configurations 
that overlaps with the Llama ones considered above, i.e. 
X = (±0.8, ±0.6, ±0.2). One sees that the result is com¬ 
patible with the previous one, though improved in the 
following aspects: (i) the NR curves extend from larger 
values of j up to the final state corresponding to the mass 
and angular momentum of the final black hole. The fi¬ 
nal state, as read from the metadata.txt file in the SXS 
catalog, is also indicated by the grey square marker in 
the plots; (ii) due to the improved accuracy of NR data, 
the differences between EOB and NR merger states (in¬ 
dicated by markers in the plots) are visibly smaller (and 
more homogeneous) than with Llama data. Compare in 
particular the case ±0.2. Generally speaking, the EOB 
and NR curves are seen to progressively separate as the 
spin becomes negative and large. Still, the final merger 
states are closer when the spin is negative than when 
the spin is positive. Such differences are typically larger 
than the numerical error, when available (shaded region 
in the plot): this suggests that the performance of the 
EOB model should be improved close to merger when 
the BHs are spinning. 

The insets in the panels of Fig. [14] are close-ups of the 
initial junk-radiation-transient parts. Note that the mag¬ 


nitude of the junk radiation transient (and thus of the 
shift (A j°, A E^) that has to be applied) is very different 
depending on the configuration. The cases x = ±0.2 are 
particularly remarkable. For them, rather large vectorial 
shifts need to be applied to reconcile the NR curve with 
both the EOB at low frequencies and the final state. 

The difference between Llama and SXS merger states 
is rather small in absolute value (~ 10 —4 ), and compati¬ 
ble with the expected level of conservation of energy and 
angular momentum in the Llama code (see Table 2 of 
Ref. [51] j. Note that the energy conservation in SXS data 
is on average more accurate by two orders of magnitude 
(see Table II of Ref. [22]). The compatibility among NR 
merger data is highlighted in Fig. [15] the figure plots the 
functions E™ s (x) and j mrg (x) (including also extremal 
spin values, see below): Llama and SXS data are hardly 
distinguishable on this scale. The SXS NR merger states 
E™ rg (x) and j mrg (x) shown in the picture can be accu¬ 
rately fitted with quadratic polynomials in x, 

Ksxsix) = -0.035546x 2 - 0.070311x - 0.119444, 

(23) 

jgx i(x) = -0.180948X 2 + 1.068176X + 2.805593. 

(24) 

The figure also summarizes the performance of the EOB 
model at merger (the corresponding points are indicated 
by red, empty, squares), highlighting that the EOB/NR 
agreement between merger states, which is best when 
X = 0, slightly worsens when \ < 0 and worsens more 
markably when \ > 0, and especially when x -+ 1. Such 
a disagreement calls for improvements in the analytical 
EOB model that will be undertaken in future studies. 
We recall however that, despite such discrepancies in the 
energetics, the phasing provided by the EOB model dis¬ 
cussed here is in excellent agreement with NR data, as 
illustrated by the dephasing and unfaithfulness data of 
Table U 


2. Extremal spins and comparison with the results of 
Taracchini et al. m 

Let us finally analyze in detail some (equal-mass, 
equal-spins) quasi-extremally spinning binaries that are 
present in the SXS catalog, x = —0.95, \ = +0.98 and 
X = +0.994, the highest spin value simulated so far. We 
recall that before the SXS catalog became public, data 
with x = —0.95 and x = +0.98 were used by Tarac¬ 
chini et al. m to test the phasing and energetics of their 
EOBNR model. We note that the EOB model of Ref. [TT] 
differs in many aspects from the one we are presenting 
here. Since we have presented a new computation of 
Eb(j) taking advantage of an improved understanding 
of various subtleties, we think it is pedagogically use¬ 
ful to compare and contrast our new NR-based results 
with both the EOB and the NR curves of Ref. JTT]. The 
latter data were kindly given to us by A. Taracchini and 
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FIG. 14. Energy versus angular momentum curves for equal-mass binaries with spin x either aligned or anti-aligned to the 
orbital angular momentum as obtained from SXS dataset. This cross checks the above comparison with Llama data. 


A. Buonanno. The result of our comparison is illustrated 
in Fig.[l6] Each top panel of Fig.[l6]contrasts four curves: 

(i) our newly computed NR curve, using a suitable vecto¬ 
rial shift, (black); (ii) our new EOB curve (red, dashed); 
(iii) the EOB curve of Ref. [IT] (light blue); (iv) the NR 
curve of Ref. HD (dot-dashed, magenta). In addition, in 
the bottom panel of Fig. [16] we show the corresponding 
differences together with the NR numerical uncertain¬ 
ties, depicted as a shaded grey region. Such NR uncer¬ 
tainties were computed by the following procedure: (i) 
for the two highest resolutions available, we determined 
the corresponding E b (j) curves and the related vectorial 
shifts (Ajo, AEfi); (ii) we took the difference between 
the so computed highest and second highest resolution 
data, AE b NRNR ; (iii) the shaded region corresponds to 
±A£^ rRNR . The uncertainty we obtain in this way looks 
rather small: it is ~ ICR 6 for large values of j, to grow 
up to a mere ~ 10~ 4 towards merger. 

A careful inspection of this figure tells us that: 

(i) The energetics yielded by our EOB model is con¬ 
sistent with extremally spinning NR data, and ex¬ 
hibits the continuation of the trend we found above 


for |x| < 0.8 . For both spins, our EOB curves are 
slightly above the NR ones close to EOB merger. 

(ii) Regarding the comparison with the EOB model of 
Ref. m we note first that the two EOB’s are per¬ 
fectly consistent among themselves (and with the 
NR data) during most of the early inspiral. The 
two EOB models start differing among themselves 
(and with NR) when getting close to merger. For 
X = —0.95 the EOB model of Ref. STT] gives a bet¬ 
ter approximation to the final state than our model; 
on the contrary, for x = +0.98 our model performs 
quantitatively better near merger. Note that, in 
spite of these differences in energetics, the wave¬ 
forms delivered by both EOB models agree well 
with NR data within the NR uncertainty HU 

(iii) On the same plot, we superpose the NR E b (j) 
curves presented in Ref. [Tl] (dashed lines; red 
online). The curves exhibit significant differences 
with both EOB models already in the early inspi¬ 
ral. Moreover their end point (after merger and 
ringdown) differs from the actual final mass and 
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FIG. 15. Equal-mass, equal-spin case. Comparison between 
EOB and NR merger quantities, binding energy (top) and 
angular momentum (bottom). 


angular momentum of the BH (as extracted from 
the SXS file metadata.txt and displayed in the 
figure). By contrast we note that we determined 
the additional, needed, final shifts taking care that 
this last point (computed with high accuracy) ap¬ 
proximately ( 1CT 4 level) coincides with the final 
point of the E^{j) curve. No details were given in 
Ref. El on how their (j ) curve was computed. 

(iv) As a side remark, we note that the red and light- 
blue curves in each insets of Fig. [16] differ at the 
~ 10 -5 level. Since both curves were obtained sub¬ 
tracting the corresponding EOB to the same NR 
curve, this 10 -5 level essentially quantifies the an¬ 
alytical difference, in the early inspiral, between the 
conservative dynamics of our EOB model and that 
of SEOBNRv2. 


To conclude, Fig. 17 illustrates the EOB/NR compar¬ 
ison for the highest value of the spin available, x = 
+0.994. The numerical uncertainty, computed as out¬ 
lined above, is again rather small, as it does not ex¬ 
ceed the 10 -4 level towards merger. The performance 
of the EOB in this case is essentially analogous to the 
X = +0.98 case illustrated above, though with a slightly 
larger displacement between the EOB and NR merger 
states. 



FIG. 16. Energetics in the extremally spinning case: com¬ 
parison with the EOB and NR data of Ref. HU. Markers 
indicate the merger for each dataset (except the NR data of 
Ref. HU): NR (red); EOB model of Ref. HU (lgiht-blue); 
our EOB model (red). The shaded area gives an estimate of 
the numerical uncertainty. Both EOB models are consistent 
among themselves and with the NR curves computed in this 
paper. For \ = +0.98 our EOB model is closer to NR than 
the EOB model of El , while the reverse is true for \ = —0.95. 


3. Unequal-mass, spinning binaries 

We conclude this survey of EOB/NR energetics com¬ 
parison by also discussing explicitly a few representative 
configurations with both unequal masses and unequal 


spins. The comparison is illustrated in Fig. [T8j which col¬ 
lects six representative configurations, notably the high- 
mass-ratio ones (8, ±0.50, 0). We show here explicitly 
those cases where only one of the two black holes is spin¬ 
ning, since in these cases we have enough resolutions to 
compute error bars. The two cases (3, —0.50, —0.50) and 
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FIG. 17. EOB/NR comparison for q = 1, xi = X 2 = +0.994 
(bottom), the highest spin simulated so far. The NR uncer¬ 
tainty (shaded region) is obtained by taking the difference 
between the data at the two highest resolutions after deter¬ 
mining the corresponding vectorial shifts (Aj°, AE°). 


(3, +0.50, +0.50) yield equivalent results. The numerical 
uncertainties, computed as in the previous section, are 
represented by light-gray shaded areas: notably, these 
error bars are compatible with the extremally spinning 
cases discussed above, i.e. about 2 x 10~ 4 towards merger 
(we expect to get similar results also for the other config¬ 
urations). Generally speaking, we see here the same fea¬ 
tures we were finding in the equal-mass, equal-spin cases 
shown in Figs. [T4}[T6| and[T7l I n particular, one sees that: 
(i) when the BH spin is positive, the EOB model predicts 
values of the merger states that are slightly smaller than 
what they should be. This occurs independently of the 
mass ratio. Similarly, (ii), when the BH spin is nega¬ 
tive, the EOB predicts values that are slightly, but sig¬ 
nificantly, higher than the NR prediction. If both black 
holes are spinning, these effects remain qualitatively the 
same, though quantitatively are slightly magnified. In 
other words, the Eb(j) analysis seems to indicate that the 
current EOB model is slightly overestimating the magni¬ 
tude of the spin-orbit interaction. Note however that the 
E h (j) diagnostics is here magnifying a difference that cor¬ 
responds to a very small part (near merger) of the time- 
domain waveform. Let us only note here that this shows 
the usefulness of the Eb(j) diagnostics for understanding 
subtle physical effects in the spin-orbit interaction that 
cannot be appreciated by just inspecting the phasing. 


VII. PRELIMINARY COMPARISON OF THE 
TWO EXISTING SPINNING EOB MODELS 


Let us finally highlight some differences between the 
analytic structure of our EOB model and that of Ref. m 
by comparing some crucial quantities that enter the con¬ 
struction of the models. 


As a first example, we focus on the q = 1, nonspinning 
case. Figure 19 compares three A(r) curves: our NR 
tuned A or b(r) curve, with ag given by Eq. §; the A(r) 
of Ref. m , calibrated to the same sample of nonspinning 
NR waveform data we used here; and the Schwarzschild 
A(r). The markers on the plot indicate the location of 
the (adiabatic) “light-ring”, defined as the peak of the 
effective “photon-potential” A(r)/r 2 (corresponding to 
r = 3 M for the Schwarzschild metric). Let us first 
note a qualitative difference: our choice of resumming 
the Taylor-expanded A function with a Pade (1,5) ap- 
proximant entails that our function A or b(r) vanishes at 
r = 0. This is qualitatively different from A Tar (r), whose 
(non Pade) resummation is chosen so as to impose the 
existence of a horizon at a nonzero value of r. As a 
consequence, our radial potential is more repulsive than 
A Tar (r) when r < 2. 


From the quantitative point of view, let us note that 
our A or b(r) and A Tar (7') are rather close up to the 
lightring of A Tar (r). Their difference (up to their light¬ 
ring location, blue dot) ranges between —3 x 1CP 3 and 
approximately 7 x 10 -3 at the lightring point. Although 
such differences look small, let us emphasize that, as dis¬ 
cussed in Ref. [55], these differences are large enough to 
make the two potentials belong to different “equivalence 
classes”, i.e. to yield significantly different conservative 
dynamics. Indeed, Ref. [59] found that one needed dif¬ 
ferences smaller than the 10 -4 level to lead to indistin¬ 
guishable dynamics. 


As a second example of comparison between our EOB 
dynamics and that of Ref. m we compare in Fig. [20|the 
dimensionless effective potential H e g(r,£,x) (including 
spin-orbit and spin-spin interactions, but with p Tt = 0), 
for x = Xi = X 2 = +0.98 and two representative val¬ 
ues of the dimensionless orbital angular momentum £. 
For the value £ = 2.40, both effective potentials exhibit 
a local minimum, corresponding to a stable circular or¬ 
bit. On the other hand, for l = 1.81 (which corresponds 
to the merger location of the blue curve in Fig. 19) the 
model of Taraccliini et al. no longer allows for a stable 
circular orbit, while ours still does. In other words, our 
dynamics is still inspiralling, while that of Ref. m is al¬ 
ready “plunging”. Note that for x = Xi = X 2 = +0.98 
we have seen above that our energetics is closer to the 
NR one. It should, however, be kept in mind that we are 
here comparing two conservative dynamics while the real 
EOB dynamics is also (especially near merger) strongly 
modified by radiation reaction effects. 
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FIG. 18. Energetics comparison for a selected sample of unequal-mass, unequal-spin binaries. The EOB/NR compatibility is 
essentially comparable to the equal-mass results shown in Figs. |11| |16| |17| 


VIII. CONCLUSIONS 

In this paper we have presented a new NR completion 
of the spinning, nonprecessing, EOB model introduced 
in Ref. [T]. 

Our new EOB model uses only two flexibility!^] param¬ 
eters, one for the nonspinning sector a§(z/), and one for 
the spinning sector 03 ( 01 , 02 ; v). We have calibrated the 
functions a%(y) and 03 ( 01 , 02 ; v) by comparing EOB pre¬ 
dictions to a large sample (~ 50) of NR data waveforms 
computed with two different numerical codes (Llama 
and SpEC). This calibration is done (for every step) 
through a simple one-parameter search by requiring that 
the EOB/NR time-domain phase difference (after align¬ 
ment) stays within NR uncertainty up to merger. This is 


Let us recall that the idea of using the “flexibility” of the 
EOB formalism to determine effective values of yet uncomputed 
higher-order analytical functional parameters (such as ag(i')) by 
fitting EOB predictions to NR data was introduced, and first 
implemented, in Ref. m- 


achieved in two steps: the function a%{v) is determined 
first by using nonspinning data; then the single remaining 
flexibility parameter C 3 (v; a 1 , 0 , 2 ) is determined by using 
spinning data. 

The EOB/NR comparison is done both for phasing 
and energetics (gauge-invariant relation £/(/) between 
binding energy and angular momentum). 

Our main results are: 

(i) The i = to = 2 GW phasing performance of our 
new spinning EOB model (measured both by 
the time-domain EOB/NR phase difference and 
through faithfulness computations) is at least as 
good as that of the most recent version of the spin¬ 
ning EOBNR model developed by the group of 
A. Buonanno nja im hbj. To quantitatively as¬ 
sess the quality of our unfaithfulness results, we 
display in Fig. [2l] a histogram of the maximized- 
over-mass values of the unfaithfulness]/] F max = 


' We recall that those unfaithfulnesses are only integrated over the 
NR frequency range. 
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FIG. 19. Comparison between EOB radial potentials A(r)\ 
the Schwarzschild case is contrasted with the q = 1, nonspin¬ 
ning function of Taracchini et al. El and the Pade-resummed, 
NR completed one of this paper, Eq. with a%(v) given by 
Eq.©. The marker indicates the location of the “light-ring”, 
i.e. the peak of the effective “photon-potential” A(r)/r 2 (lo¬ 
cated at r = 3 M for Schwarzschild). 



FIG. 20. Comparison between effective Hamiltonians 
Hes(r,l,X ) obtained by using our EOB spinning model and 
the one of Ref. in for x = +0.98 and for two values of the di¬ 
mensionless orbital angular momentum t. The value i = 1.81 
corresponds to the merger location of the blue curve in Fig. 1 19| 
For i = 1.81 our effective potential is such that the system is 
still inspiralling. 


maxjvf (F(M)), over the entire sample of SXS wave¬ 
forms we use. Note that the histogram is skewed 
towards very small values of F maxi with a me¬ 
dian value equal to 0.00056105 (see vertical line on 
the figure). It is interesting to note that F max < 
0.003 for all data except for the case (q,x I 1 X 2 ) = 



FIG. 21. Maximum value of the EOBNR unfaithfulness F 
(listed in Table |I| for all SXS dataset considered. The vertical 
line indicates the median of the histogram, corresponding to 
0.00056105. 


(1, —0.50, 0) and the near-extremal spin q = 1 data 
sets with x = Xi = X 2 = (+0.98,+0.99,+0.994). 
For the latter datasets, not that \ = +0.99 and 
X = +0.994) both yield F max ss 0.005, while 
X = 0.98 yields F max « 0.0042. 

When comparing the mass dependence of F in 
our Fig. [9] with the corresponding quantity dis¬ 
played in Fig. 1 of Ref. HU one notices the fol¬ 
lowing facts: (i) our F(M) is, for all values ex¬ 
cept a few, maximal around M = 20M 0 and then 
monotonically decreases as M increases. By con¬ 
trast, the mass dependence of the F{M) of Ref. [TT| 
is quite different: either it increases monotonically 
from 20 Mq to 2OOM 0 or it has a local maximum. 
We attribute the mass behavior of our F(M) to 
the new, more accurate attachment of the ring- 
down part [l4j which we have implemented in our 
model; (ii) In Fig. 1 of HU we spot two outliers with 
F max ~ 0.01: one is SXS:BBH:0065, with q = 8 
(Xi,X 2 ) = (+0-5,0); the other outlier of Ref. [TT] is 
SXS:BBH:0152, with q = 1, Xi = X' 2 _= +0.6. We 
stress that our corresponding values F(M) remain 
below 10 —3 , as displayed in Fig. [ 9 } 

(ii) We have presented the first systematic, accurate 
computation of E^j) curves from SXS and Llama 
data of nonprecessing, spinning (or nonspinning) 
configurations. Our results significantly improve 
previous attempts to do so and give us a reliable 
target for EOB calibration, that complements the 
usual phasing analysis. 
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For the nonspinning sector of the model, the in¬ 
spection of the NR Et(j) curves led us to decide 
that the choice T rt = 0 (which was already used in 
early EOB works [35112S]) leads to a good EOBNR 
agreement between the corresponding energetics, at 
the 10 -4 level up to merger. We hope to further 
study (both analytically and numerically) the role 
of the radial part of radiation reaction, and its in¬ 
fluence on Schott terms. The results presented in 
Figs. [4] and [5] have shown the potential of SXS 
data for very accurate investigations of energetics 
in coalescing black hole binaries. 

The Eb(j) NR curves were computed using two dif¬ 
ferent sets of numerical data, one by the Llama 
code and the other by the SpEC code. We solved 
several subtle issues related to the computation 
of Eb(j) from SXS data. In particular, for the 
cases of extremal spins Xi = Xi — —0.95 and 
Xi = Xi = +0.98, we showed that the correspond¬ 
ing NR curves presented in Ref. m were inconsis¬ 
tent both with the NR final state (as provided in 
the SXS files itself) and with early inspiral EOB 
predictions. This allowed us to meaningfully com¬ 
pare, for the first time, extremal-spin NR data both 
to our EOB model and to the one of Ref. mi. 
The main outcomes of this comparison are: (i) for 
X = +0.98, the energetics predicted by our model 
model is closer, than that of the model of Ref. mi. 
to the NR one, during the late inspiral and near 
merger; (ii) for x = —0.95, the opposite is true, 
with the model of Ref. mi being closer than ours 
to the NR values near merger. 

(iii) The EOB model we are presenting here is analyti¬ 
cally simpler than the one of Ref. mi, while it per¬ 
forms at a comparable accuracy (which is demon¬ 
strably better in some situations, e.g. the case 
of SXS:BBH:0152 data). The structural simplic¬ 
ity of our spinning EOB model is such that it is 
easy to recalibrate it if the needs occurs, e.g. in 
the case that further linear-in-^ corrections to the 
gyro-gravitomagnetic ratios ( Gs,Gs „), as defined 
in Paper I, are introduced. Different analytical ex¬ 
pressions would imply a different determination of 
the function 03 ( 0 , a; 11 ). In addition, it is straighfor- 
ward to incorporate tidal effects in the model sim¬ 
ply by augmenting the pure orbital A function with 
the tidal potential recently determined in Ref. m 
so as to build a BNS (or BHNS) spinning EOB 
model. The performance of such a tidal, spinning, 
EOBmodel, as well as its eventual calibration, will 
have to be carefully assessed against state-of-the- 
art NR simulations of spinning NS binaries 1531 . 

(iv) The preliminary comparison given here between the 
performances of our model and the EOB model of 
Ref. mi , was purely based on information that was 
either given to us or that was extracted from the lit¬ 
erature. We also emphasized some of the structural 


differences between these two spinning EOB mod¬ 
els: (i) in Fig. 19lwe constrasted the two different 


EOB A(r) radial potentials for q = 1 nonspinning 
black holes; and (ii) in Fig. 20 we contrasted the 
two effective Hamiltonians H e s(r) for q = 1 and 
X = +0.98. The differences are nonnegligible and 
are probably partly responsible of the different per¬ 
formances of the two EOB models versus the NR 
Eb(j) curve for x = +0.98. We hope that it will be 
possible to perform soon a more detailed compari¬ 
son between the two spinning EOB models, notably 
by directly comparing waveforms. We think that a 
synergetic effort dedicated to transfering the best 
features of each EOB model into a new model for 
spinning binaries will be crucial for the forthcoming 
gravitational wave astronomy. 
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Appendix: Computing using SXS waveforms 


In this Appendix we discuss in detail the procedure we 
adopt to compute the function Eb(j) using SXS NR data. 
First, in Sec. [l] we contrast SXS with Llama data and 
highlight some (minor) unsatisfactory features of the SXS 
waveforms. The actual Eb(j ) computation is detailed in 
Sec. m below. 
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1. Unphysical effects in the q = 1, SXS multipolar 
waveforms 


For an equal-mass binary, with equal (aligned or anti¬ 
aligned) spins, the fact that the system is symmetric by 
exchange of the two black holes ensures that all multi¬ 
poles with to = odd have to vanish. We verified that this 
is essentially the case for all Llama simulations at our dis¬ 
posal: in the time plots of these multipoles, only uncorre¬ 
lated noise is found. On the contrary, a careful inspection 
of SXS data led us to the discovery that the m = odd 
multipoles have amplitudes, which are small but which 
are not just uncorrelated noise. Rather they qualitatively 
show some sort of chirping structure present in the ac¬ 
tual physical modes with to = even. This is illustrated in 
Fig. [22] for the case q = 1, = Xi = +0.6, correspond¬ 

ing to dataset SXS:BBH:0152 of the SXS catalog. The 
figure contrasts the SpEC modes (top panel) with the cor¬ 
responding Llama modes (bottom panel). One immedi¬ 
ately sees the striking qualitative difference between the 
two independent simulations of the same physical system. 
Though, the amplitude of these modes is so small that it 
does not have any practical influence on the Ef,(j) com- 
putatiorj^] it is interesting to briefly point out the impact 
of two physical effects on shaping the multipolar struc¬ 
ture of the modes: (i) the motion of the center of mass 
(CoM) of the system and (ii) the (tiny) asymmetry of the 
(nominally symmetric) initial binary configuration. 

The effect of the motion of the CoM due to residual lin¬ 
ear momentum in the initial data has been investigated 
in some very recent work [5U [55]. Notably, Ref. [55] 
introduced a method to eliminate such residual momen¬ 
tum from the initial data, so that future SXS waveforms 
will be amended of such systematic uncertainty. How¬ 
ever, as emphasized by Boyle |54j , the waveforms cur¬ 
rently present in the SXS catalog are not expressed in 
the CoM frame (as are the EOB ones) and this intro¬ 
duces systematic uncertainties that should be removed. 
As a matter of fact, Ref. )54j (that appeared while this 
paper was under evaluation) proposed a way to remove 
the effect of the CoM drift from the waveforms. We asked 
Micheal Boyle to kindly apply his removal procedure to 
the SXS:BBH:0152 configuration we are discussing here. 
His result is shown (with his permission) in Fig. 23 One 


sees that the secular growth of the mode amplitudes is 
visibly reduced (though it is not completely removed, es¬ 
pecially towards merger). This is particularly visible for 
the (2,1) mode, where the trend visible in Fig. 22 is re¬ 


than before, and visually more consistent with the out¬ 
come of the Llama simulation. Improvements in the de¬ 
termination of the CoM frame of a NR simulation may 
flatten the to = odd modes even further. 

Let us now turn to estimating the effect on the to = 
odd multipoles due to the tiny asymmetry present in 
the NR initial data. Inspecting the file metadata.txt of 
dataset SXS:BBH:0152, one finds that the initial massed 
and initial spins are actually slightly different; i.e., one 
has toi = 0.499999998745987, to 2 = 0.499999997467155 
and xi = 0.599963311735964, y 2 = 0.599963312827750. 
These values of the initial masses yield v = 0.25 at ma¬ 
chine precision, though Xi — X 2 = {mi — to 2 ) /M = 
1.278831984752316 x 10~ 9 . We then estimated the im¬ 
pact of these tiny asymmetries by running an EOB simu¬ 
lation with precisely these spin parameters and this value 
of Xi — X 2 entering the subdominant odd-parity modes. 
The multipolar amplitudes of the EOB waveform com¬ 
puted with these choices are exhibited in Fig. [24} one 
sees that during inspiral, where the analytical (EOB, re¬ 
summed) waveform is highly reliable, the amplitude of 
subdominant modes is from 1 (for the ip 2 \ mode) to ap¬ 
proximately 2 (for the +43 mode) orders of magnitude 
smaller than the corresponding SXS multipoles in the 
CoM frame. In Fig. 24 we also added the SXS (2,1) 


placed by random-like oscillations. Globally, the ampli¬ 
tude of the modes is rather constant, with less structure 


mode, in the CoM frame, to ease the comparison. This 
finding suggests that the magnitude of the subdominant 
multipolar amplitudes we see in Fig. [23] is only partially 
due to the asymmetry in the initial data and that other 
effects are present. In any case, since the analysis we 
did here confirms that the effects of the asymmetries are 
practically negligible, we are entitled to ignore them and 
consider exactly symmetric initial data for the EOB evo¬ 
lutions. 

Finally, let us note in passing that ~ |'033° B | 

up to merger, likewise the corresponding SXS multi¬ 
poles in the CoM frame in Fig. JH This result con¬ 
trasts with the top panel of Fig. 22 where the (2,1) and 
(3, 3) modes were not coinciding, and gives further ev¬ 
idence of the need of removing the effects of the CoM 
motion motion from the SXS waveforms. Although the 
numerical effects we are pointing out here are small and 
practically negligible for the equal-mass configurations, 
they may not be as small in other situations (e.g., as 
pointed out in Ref. [51] for the SXS:BBH:0004 dataset, 
( 91 X 11 X 2 ) = (1,—0.50, 0). In principle, we think that 
a careful analysis of the (hopefully very small) impact 
of the CoM drifts on all the EOB/NR comparisons dis¬ 
cussed in this paper, both for phasing and energetics, will 
be needed in the future. 


We have analyzed what happens if the m = odd modes are put to 
zero. We find that the resulting small difference (mainly coming 
from the m = odd junk radiation) can be absorbed by modifying 
the vectorial shift (Aj°, AE®) mentioned in the main text and 
discussed below. In practice, all our E(j) curves are computed 
using the full multipolar information at our disposal including 
multipoles up to i = m = 8. 


2. Computation of Eb(j) from SXS waveform data 

Let us now turn to explaining how Eb(j) curves are 
computed in practice. First of all, we need to iden¬ 
tify the initial state of the system, i.e. the initial 
(^adm^adm)- For each simulation, this information 
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FIG. 22. Analysis of the subdominant multipoles of q = 1 
SXS waveforms (top panel). The figure compares the ampli¬ 
tudes of various subdominant multipoles to the l = m = 2 
one. The multipolar amplitudes with m = odd grow in time 
likewise the dominant one. This effect (that is found in all 
q = 1 data) is (partly) related to the motion of the center 
of mass and to residual tiny asymmetries in the intial data, 
as discussed in the text. It is, however, absent in the corre¬ 
sponding Llama waveform data (bottom panel). 


is provided in the file metadata.txt that is found in 
any of the Lev* directories that can be downloaded 
at http://www.black-holes.org/waveforms/ for any 
given configuration. Then, the next step is to take all 
multipoles up to £ = 8 and use them to compute the 
fluxes of energy and angular momentum as functions of 
time. In inspecting them, one quickly realizes that all 
simulations basically split into two equivalence classes 
depending on the properties of the initial transient in 



FIG. 23. The same SXS multipolar waveform of Fig. |22| after 
M. Boyle removed the drift of the CoM using the procedure 
of Ref. [5j|] . The subdominant multipolar amplitudes look 
flatter than in Fig. [22] and qualitatively closer to the Llama 
ones. 



FIG. 24. Investigating the effect of tiny asymmetry in the 
initial data at the level of EOB waveform. The (2,1) mode 
of Fig. [23] is also included for visual comparison. The asym¬ 
metry in the initial data is not sufficient to fully explain the 
magnitude of the subdominant modes in the data of Fig. |23| 


the fluxes driven by the junk radiation: in class (A), one 
sees a huge burst of energy and angular momentum radi¬ 
ated during approximately the first 150 M of evolution. 
This behavior is displayed in the top panel of Fig. [25j 
that refers to the \i = Xi = +0.98 case. Once the total 
losses at time t are obtained (by integrating the fluxes 
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FIG. 25. Computation of Eb(j) for \i = X 2 = +0.98. The 
huge burst of junk radiation (top panel) introduces a huge, un¬ 
physical, vectorial shift of the raw curve (grey, dashed), whose 
end point is strongly displaced from the actual final state 
(jf, Mf ) (grey square on the plot) given in the metadata.txt 
file of the SXS catalog. The solid curve is obtained by ap¬ 
plying a shift vector (A j°, A E°) to the raw curve so that the 
final point of the shifted Eb(j) coincides with ( jf,Mf ). The 
filled circles mark the merger location. See text for further 
details. 


FIG. 26. Different phenomenology for = X 2 = +0.60. 
The initial burst of radiation is much more moderate than 
the x = +0.98 case so that the direct, raw computation of 
Eb{j) looks already visually consistent with the final state 
and the EOB curve. Additional fine tuning of (Aj°,AE°) 
is still allowed and it is then performed to further minimize 
the EOB-NR difference (within the NR uncertainty) for large 
values of j. See text for details. 


up to time t), the presence of such a big burst of ra¬ 
diation results in a huge shift in the (j, E^) plane away 
from (Jadmi-^adm) 1 This, raw, curve is represented as 
a gray, dashed line in Fig. [25] The location of the merger 
is indicated by a filled, grey, circle. The unphysical na¬ 
ture of the initial, junk-related, shift is indicated by the 
fact that the last point of the Ef,(j) curve, which should 
coincide, to a good approximation, with the mass and 
angular momentum of the final black hole ( jf,Mf ), sig¬ 
nificantly differs from them. 


The value of ( jf,Mf) is among the data given 
in the SXS catalog, and is notably included in the 
metadata.txt file. In Fig. [25} the actual ( jf,Mf) is 
represented by a filled, grey, square on the plot, while 
the end point of the raw Eb(j) curve (indicated by an 
empty square) is noticeably displaced from ( jf,Mf ) due 
to the effect of the junk radiation. Such displacement 
should be corrected, as we shall discuss below. In the 
second class of data, that we call (B), such an unphysical 
initial burst of GW radiation is much smaller and the 
raw computation of the E b (j ) curve gives a result that, 
without any additional shift, is visually consistent with 
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FIG. 27. Choosing A j° by inspecting the EOB/NR dif¬ 
ferences between dimensionless orbital frequencies fl for two 
representative datasets. The bottom subpanel of each con¬ 
figuration illustrates how Af2 EOBNR is flattened (and the 
eccentricity-driven oscillation averages zero) on a large j- 
interval once a good choice of A j° is made. 


the known final NR state. This is illustrated for the case 
Xi = X 2 = +0.60 in Fig. [26j Still (as highlighted in the 
inset of the figure) the NR curve, when zooming on the 
difference in the early part (see inset), exhibits a signif¬ 
icant early deviation from the EOB, contrary to what 
was happening for the xi = X 2 = +0.60 case with Llama 
data, a simulation starting at a smaller value of j. 


TABLE III. The shift vector (Aj°, AE°) that is applied to 
the raw Eb(j) curves computed for an illustrative subset of 
the configuration listed in Table [T] The numbers listed refer 
to curves where the final state (j f , Mf ) was already imposed 
on the raw data, except for those datasets marked with an 
asterisk, *. 


Name 

9 

Xi 

X2 

A E° b 

O 

<1 

SXS:BBH:0066 

1 

0 

0 

—5.583e-03 

—2.13e-02 

SXS:BBH:0169 

2 

0 

0 

+1.224e-05 

—5.0e-04 

SXS:BBH:0030 

3 

0 

0 

—1.669e-03 

—9.0e-03 

SXS:BBH:0167 

4 

0 

0 

+1.917e-04 

+7.0e-04 

SXS:BBH:0056 

5 

0 

0 

—2.571e-03 

— 1.25e-02 

SXS:BBH:0166 

6 

0 

0 

+1.942e-04 

+7.0e-04 

SXS:BBH:0063 

8 

0 

0 

—7.7709e-04 

—4.0e-03 

SXS:BBH:0185 

9.989 

0 

0 

+4.053e-04 

—2.4e-03 

SXS:BBH:0156 

1 

-0.95 

-0.95 

—9.496e-05 

—6.0e-03 

SXS:BBH:0154* 

1 

-0.80 

-0.80 

—3.018e-04 

—0.9e-03 

SXS:BBH:0151* 

1 

-0.60 

-0.60 

—6.507e-05 

—5.7e-04 

SXS:BBH:0149 

1 

-0.20 

-0.20 

—4.460e-04 

—3.8e-03 

SXS:BBH:0150 

1 

+0.20 

+0.20 

—4.200e-04 

—2.8e-03 

SXS:BBH:0152* 

1 

+0.60 

+0.60 

—7.393e-05 

—0.7e-03 

SXS:BBH:0155* 

1 

+0.80 

+0.80 

—2.536e-04 

—0.5e-03 

SXS:BBH:0172 

1 

+0.98 

+0.98 

—3.640e-03 

— 1.05e-02 

SXS:BBH:0178 

1 

+0.994 +0.994 

—2.657e-03 

—8.8e-03 

SXS:BBH:0162 

2 

+0.60 

0 

—5.240e-04 

—2.3e-03 

SXS:BBH:0036 

3 

-0.50 

0 

—2.434e-03 

—5.5e-03 

SXS:BBH:0031 

3 

+0.50 

0 

—2.326e-03 

—7.1e-03 

SXS:BBH:0064 

8 

-0.50 

0 

—4.604e-03 

—5.5e-03 

SXS:BBH:0065 

8 

+0.50 

0 

—4.792e-03 

—6.8e-03 


Finally, in both cases, (A) or (B), we need to apply a 
shift vector (Aj°, A E°), which will be larger in case (A). 
In a first step (which is mostly useful for case (A)) we ap¬ 
ply a first shift (Ajf, AEjj) allowing the final point of the 
Et(j) curve to essentially (within < 10 -4 ) coincide with 
(jf. Mf). In a second step, we refine the determination of 
the vectorial shift by performing a second, smaller, shift, 
(Aj°,AF“). 

The second shift vector, simply denoted (Aj°, AE°), 
is determined as follows. First of all we look at the 
dimensionless “orbital frequency” fl = Mfl, that is de¬ 
fined (both in the EOB model and NR) as fl = dE^/dj. 
This quantity is, by construction, independent of the con¬ 
stant A E$, but depends on A j°. We then compute the 
difference Afi EOBNR (j) = f) EOB (j) — 0 NR (j) and plot 
it versus j. The shift A j° is chosen so that Afl(j) is 
“flattened” as much as possible on the largest possible 
j-interval. In practice, we choose A j° such that the oscil¬ 
lating (because of residual, tiny, eccentricity effects in the 
NR waveform) A fl(j) averages to zero. Figure 
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shows 

the oscillating Afl(j) for two, fiducial, configurations cho¬ 
sen randomly in the data sample: (q, X 11 X 2 ) = (2,0.6, 0) 
or (q,Xi>X2) = (3,—0.50,0). In both cases, the raw 
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AQ (blue line) shows a trend, that can be flattened 
by just choosing properly the constant A jo (red curve). 
Once this is done, we inspect the difference AI5 BOBNR = 
E BOB — E^ R and we find it remains practically constant 
on a large j interval. We then choose A E® so that the 
AE BOBNR is found to oscillate around zero. This algo¬ 
rithm for determining the shift vector (Aj 0 , A E®) can be 
applied to all SXS spin-aligned configurations at our dis¬ 
posal. Note that the shifts that are so determined are, 
in most cases, tiny, but relevant on this scale. In some 
cases (see e.g. case (q,xI1X2) = (8, — 0.5,0) in Fig. 14) 
the combination of the two shifts leads to an end point 
for the shifted E b {j) NR curve that is close to, but vis¬ 
ibly displaced with respect to the (j/, Mf ) provided in 
the metadata.txt. Note however that in this case the 


total angular momentum changes sign during the inspiral 
which may be connected to subtleties in the computation 
of E b (j). We leave to future work a deeper investigation 
of this and similar cases. 

For completeness, we list in Table |III| the shift vec¬ 
tors (Aj°, AFj)) that are used to improve the compati¬ 
bility between (Aj°,AE®) the NR E b (j) and the EOB 
one. The numbers listed in the table always indicate 
the crucial “second-shift” vector mentioned above (i.e., 
are applied to curves already shifted by (Ajf ,AE]i) so 
as the final point of E b (j) approximately coincides with 
( jf,Mf )), except for those datasets marked by an as¬ 
terisk, where (Aj°,AE b ) is applied directly to the raw 
data. 
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